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analytically according to the macroscopic wave functions of the condensates, 
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I. INTRODUCTION 

Since the realization of Bose-Einstein condensates (BECs) in dilute atomic gases, a num- 
ber of interesting experiments have been conducted to investigate multispecies Bose gases, 
in which two or more states of condensates exist together in a magnetic or optical trap [1,2]. 
Recently, vortex states have been obtained in a two-species Bose gas [3]. Progress in the 
experiments exploring dilute mixtures of quantum gases has stimulated intensive research 
on the properties of mixed Bose gases at zero temperature [4-6] and finite temperature as 
well [7]. 

The BEC trapped in an optical lattice exhibites a novel feature, namely the quantum 
phase transition between a Mott-insulator and a superfluid [8]. Such quantum phase tran- 
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sition has attracted considerable attention in recent years. As a matter of fact, atomic gas 
of bosons in BEC subjected to a lattice potential which is turned on smoothly can be kept 
in the superfluid phase as long as the atom-atom interactions are small comparing with the 
tunnel coupling. In this regime, the kinetic energy is dominant in the total energy of the 
boson system. With an increase of the potential depth of the optical lattice, it is getting 
more and more difficult for bosons to tunnel from one site to the other, and finally the 
system attends an insulator phase above a critical value of the potential depth. In this case, 
the phase coherence is absent and the number of boson atoms in each lattice site becomes 
the same. The system possesses a Mott-insulator behavior. Various approaches have been 
proposed to understand theoretically the quantum phase transition and to determine the 
phase diagram as a function of BEC parameters [9-14]. 

Motivated by both the experimental and theoretical progress, we in the present paper 
study the phase diagram for superfluid and insulator phases of two-species BECs in a one- 
dimensional (ID) optical lattice and the property of persistent current as well. The paper 
is organized as follows. In Sec. II, the exact macroscopic wave functions of the condensates 
which are not in the tight-binding regime are constructed by solving the coupled nonlinear 
Schrodinger equations. In Sec. Ill, the phase diagram is determined analytically according 
to the macroscopic wave functions of the condensates, i.e., the order parameters. Finally, 
we summarize our results in Sec. IV. 



In this ID geometry, the confinement along the radial direction is so tight that the 
trap frequency u along the radial direction is much greater than the mean-field interaction 
energy. At low temperatures, the dynamics of the atoms in the radial direction is essentially 
"frozen," with all the atoms occupying the ground state of the harmonic trap with the wave 
function that 



II. THE EXACT MACROSCOPIC WAVE FUNCTIONS 
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Here the extension of the wave function in the radial direction is given by the length scale 
lo = \Jh/mujQ of a harmonic oscillator, where m is the mass of the atoms. 

In the mean-field regime, Iq is much greater than the radius of the interatomic potential 
R e . The scattering of atoms in this effective ID system is thus still a process of three- 
dimension. According to ref. [15], the effective coupling constant in this ID system is 
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where a is the s-wave scattering length and C is a numerical constant of the order unit. The 
term Ca in Eq. (2) is negligible for l Q >> R e . In this limit, the expression for g 1D is the 
same as that obtained by averaging over the radial wave function (1), 
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We use this expression in the rest of this paper. 

We consider the two-species BECs in a ID periodic potential. The energy functional is 
seen to be 
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where ipi, and U = ^Jh/miU are the macroscopic wave functions of the condensates, the 
mass, and the harmonic-oscillator lengths in the radial direction of the ith species (i = 1,2), 
respectively, ai, and a i2 denote the s-wave scattering lengths between same-species and 
interspecies collisions. Vi(x) are the periodic potentials, 

Vi(x) = V ,isn 2 (k L x, k) , (5) 

with Vo,i denoting the magnitude of potentials, where ki = 2tt/\ is the wavevector of the 
laser light and A is the wavelength, corresponding to a lattice period d = A/2, sn^x, k) is 
the Jacobian elliptic sine function with modulus k (0 < k < 1). In the limit k — 0, the Jaco- 
bian elliptic sine reduces to the sinusoid function and thus V(x) possesses a standard form 



of the standing light wave. For values of k < 0.9, the potential is virtually indistinguishable 
from a standing light wave. Finally, for k — > 1, V(x) becomes an array of well-separated 
hyperbolic secant potential barriers or wells. 

The governing equations of the trapped BECs are obtained in terms of the variational 
procedure [16], 

which leads to the coupled nonlinear Schrodinger equations 
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For the case of weakly coupled condensates in an optical lattice [17], the wave function ■?/> 
can be decomposed as a sum of wave functions localized in each well of the periodic potential 
(tight-binding approximation) with the assumption relying on the fact that the height of the 
interwell barrier is much higher than the chemical potential. We, however, do not restrict 
ourselves to the low-energy case and look for the global condensate wave functions of exci- 
tations: ipi(x, t) = 4>i{x) exp (—ifct/h), where //j (i = 1, 2) are the chemical potentials. Thus 
the spatial wave functions satisfy the stationary coupled nonlinear Schrodinger equations 
that 
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With the general form of spatial wave functions (f>i(x) written as [18] 4>i{x) = 
Ti(x) exp [iipi (x)], Eq. (8) can be separated as real and imaginary parts. We then inte- 
grate once for the imaginary part and obtain the first-order differential equations for the 
phases (pi(x), 

ViV) = 0) 
where parameters ctj (i — 1,2) are constants of integration to be determined. Substituting 
Eq. (9) into the real part obtained from Eq. (8) and integrating again, we find 
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where fa (i = 1,2) denote additional constants of integration. 
We then construct the solutions as 

rj(x) = A iS n 2 (k L x, k) + B h (11) 

where the constants Bi (i = 1,2) determine the mean amplitudes and act as the dc offsets 
for the numbers of the condensed atoms [19], and the parameters Ai (i — 1,2) are to be 
determined. 

Substituting Eq. (11) into Eq. (10) and using identities of Jacobian elliptic functions, 
we obtain eight equations for the parameters ctj, fa, fii, and A^. Eliminating fa, we find 
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For /c = 0, sn(/c L x, 0) = sin(/c L x), the solutions reduce to 
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The constants A and B are related by restrictions such that B\ > —A\ for A\ < 
B 1 > for A° > 0; B 2 > -A° 2 for A\ < and 5 2 > for A° > 0. 
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III. THE PHASE DIAGRAM 



The average particle number densities for the two species are obtained as 
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where x' = /c^x and L = hd denotes the length of the optical lattice with h = 1, 2, 3, ... This 
leads to 

AO 
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Then the macroscopic wave functions of the condensates Eq. (15) can exist only when 
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The condensate atom currents can be evaluated from the usual definition, j = 
(h / m)Im[ip* (dip / dx)] [20], with the exact wave functions Eq. (15) which are seen to be 
travelling matter waves. The result is 
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which are independent of space-time variables and therefore persistent currents. We may 
demand that the wave functions Eq. (15) satisfy the periodic boundary condition ipi(x,t) = 
tpi(x + L,t) which is naturally fulfilled as the total length L is an integer times the lattice 
constant d (L = hd, h = 1, 2, 3, ...). These periodic solutions in ID space with spatial period 
L are equivalent to the solutions in a ring of circumstance L. The persistent currents then 
can be viewed as in the optical lattice ring. 

It is found from Eq. (22) that the persistent currents are valid only for conditions in 
which the number density of atoms is greater than critical values, n« > |A°|/2, i.e., when the 



macroscopic wave functions of the condensates exist. These persistent currents are similar 
to the ID Frohlich superconductivity induced by the traveling lattice wave [21] and can 
be controlled by adjusting the barriers height of the periodic potentials and parameters 
of the bosonic atoms. The currents increase with the decrease of \A®\ and approach the 
asymptotic maximum values j iimax = ±/mjA; £ /mj when \A®\ become vanishingly small. With 
the recent progress made on confinement of atoms in the light-induced [22,23] as well as the 
magnetic-field-induced [24,25] atom waveguides [26], the persistent currents may be observed 
experimentally in the future. 

The energy spectrum for the two species is obtained as 

n\ = E RA + ^y- + 2hu (a 1 n 1 + a 12 n 2 ), 

/il = E R>2 + + 2huj Q (a 2 n 2 + a 12 ni}, (23) 

where E Rji = h 2 k 2 L /2rrii are the recoil energy of an atom absorbing one of the lattice phonons 
[17]. 

In this paper, we consider the two species both with repulsive interactions, namely, 
di > 0, a 2 > 0, ai2 > 0. Then the macroscopic wave functions of the condensates can exist 
only when 

> E R , + % + (24) 



2 2 \a\a 2 — a 2 



121 



+ ai2 la^Vb,! - aiVo,2| 



2 \a\a 2 — a 2 



121 



H 2 > Er,2 + — + 
4 



^0,2 , a 2 |ai 2 Vo,i - aiVb,2| 



2 ' 2 \a\a 2 — af 2 | 

a 12 1^12^0,2 — a 2^,l| 



2 \a\a 2 — a 2 2 \ 

The macroscopic wave functions of the condensates ipi(x, t) are complex functions defined 
as the expectation value of the boson field operators: ipi(x,t) = (^/i(x,t)^, which have the 
meaning of order parameters and characterize the off-diagonal long-range behavior of the 
one-particle density matrix pi(x ,x,t) = (4ff(x',t)4fi(x,t)^ [16]. So the condensates can be 
kept in the superflud phase only when the macroscopic wave functions exist. Otherwise, the 



phase coherence and the currents vanish and therefore the condensates are in the insulator 
phase. Then we obtain four cases of phases for two-species BECs in a ID optical lattice as 
follows. 

Case 1. The two species are both in the superfluid phase, namely 

. a 12^0,2 — a 2^0,l , oc x 

ni > -7 ; — v 25 
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which is labeled as SS in the phase diagram, Fig. 1. 

Case 2. Species 1 is in the superfluid phase while species 2 is in the insulator phase, 

™i ^ u ( 26 ) 



n 2 < 



|ai 2 Vo,i - 01^0,2 1 



4hu \aia 2 — a\ 2 \ ' 
labeled as SI. 

Case 3. Species 1 is in the insulator phase and species 2 is in the superfluid phase, 
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labeled as IS. 

Case 4- The two species are both in the insulator phase, 
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labeled as II. 

The quantum phases of the condensates can be determined by all parameters of BECs 
and optical lattice as shown above. In Fig. 1, we show the phase diagram with the various 
same-species s-wave scattering lengths and equal particle number density n\ = n 2 = n and 
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the magnitudes of potentials (V ,i = V 0t2 = V ) for simplicity. Thus the conditions for the 
two species in the superfluid phase are given by 

and 



CtlCt2 
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respectively. 

Component separation in two-species BECs has been predicted by means of mean-field 
theory [4-7] and observed in experiments [1,2] when the relation of the scattering lengths 
that a 12 > ^Ja\(i2 is fulfilled. From the above conditions Eqs. (29), (30) and the phase 
diagram Fig. 1, we find that the larger values of a 12 favor the superfluid phase in the two- 
species mixture and the component separation according to the experimental observation 
[1,2]. Particularly when the interspecies scattering length approaches the value of the same- 
species such that a 12 = a 2 or ai, the conditions Eqs. (29), 30) result in the superfluid phase 
independent of the potential magnitude Vq. One should not be surprised by this result since 
we consider the case in which the chemical potential is always higher than the potential 
magnitude Vq seen from Eq. (24). 



IV. CONCLUSION 

In conclusion, the exact macroscopic wave functions of two-species BECs in an optical 
lattice beyond the tight-binding approximation are studied. The phase diagram is deter- 
mined analytically according to the order parameters, and persistent currents in an optical 
lattice ring are obtained explicitly in terms of the exact wave functions, which are seen to 
be traveling matter waves. 
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Figure Caption: 

FIG. 1. Phase diagrams of two-species BECs in a ID optical lattice. The magnitude of po- 
tentials Vq is in units of AnhtOQ. The interspecies scattering length is in units of ran. Dashed 
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curves: phase boundary of species 1; solid curves: phase boundary of species 2. The s-wave 
scattering lengths between the same-species are as follows: (a) a\ = 6 nm ( 87 Rb), ai = 3 nm 
( 23 Na). (b) a\ = 3 nm (hyperfine state | / = l,m/ = 1) of 23 Na), ai = 2.5 nm (hyperfine state 
| / = l,m/ = 0) of 23 Na). 
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